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Abstract 



We consider the non-equilibrium evolution of the inflaton field coupled to 
both lighter scalars and fermions. The dissipational dynamics of this field 
is studied and found to be quite different than that believed in inflationary 
models. In particular, the damping time scale for the expectation value of the 
zero momentum mode of the inflaton can be much shorter than that given 
by the single particle decay rate when the inflaton amplitudes are large, as 
in chaotic inflation scenarios. We find that the reheating temperature may 
depart considerably from the usual estimates. 
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The Universe after inflation is a cold and dark place. The inflationary expansion phase 
has driven the matter and radiation energy densities down to almost zero. If we are to 
recover the Universe of the standard hot big bang model, some source of energy (or perhaps 
more accurately, entropy) must be enlisted to rethermalize the Universe. 

Once inflation has ended, the standard picture of inflation [1],[| has it that the zero 
momentum mode of the inflaton field oscillates about the minimum of its (effective) potential. 
As it oscillates, the couplings of the inflaton to lighter particles that would give rise to single 
particle decays induce a friction-like term, which in turn gives rise to dissipational dynamics 
for the zero mode. This dissipation allows for the energy contained in the zero mode's 
oscillations to be used for the creation of light quanta which will eventually thermalize and 
allow the standard big bang picture to proceed |§. 

This scenario has been questioned recently when it was realized that the evolution of the 
zero mode of the inflaton induces large amplification of the quantum fluctuations resulting 
in copious particle production j§-§|. We call this mechanism of particle production induced 
amplification. (In the case of exactly periodic fields this reduces to parametric resonance 0). 
When the zero mode of the inflaton field has a large amplitude, this mechanism of induced 
amplification is much more efficient for particle production than single particle decay. In this 
letter we compare the time scales for damping and particle production for the mechanisms 
of induced amplification and single particle decay as well as the distribution of particles 
produced. It is important to understand the different processes, because these time scales 
determine the reheating temperature 

We consider the situation in which the inflaton is coupled to lighter scalar and fermion 
fields and track the behavior of the expectation value the inflaton zero mode, as well as 
the various modes of the lighter field. We find that the modes of the lighter field are not 
distributed thermally, but are more skewed towards modes with momentum k < M$. This 
may have implications for the thermalization phase of the reheating process. Interactions 
will tend to thermalize this spectrum to temperatures ~ M$. The case of coupling to 
fermions is qualitatively different ^from the coupling to bosons because of Pauli blocking. 

In both cases, we consider relaxation both in the linear regime, where the amplitude of 
the zero mode is not large and an amplitude expansion of the effective equation of motion 
can be used reliably, as well as non-linear relaxation, where the initial amplitude is large; 
this is typically the behavior of the inflaton in chaotic inflationary models 0. 

We use the techniques that we developed previously || to construct the effective equa- 
tions of motion for the zero mode in the presence of fluctuations due to the higher momentum 
modes. The reader is referred to this and forthcoming work || for more details. 

Throughout we will assume that the typical rates involved in the particle production 
process are much larger than the expansion rate of the Universe H; this allows us to restrict 
ourselves to the Minkowski space case. 

Linear Relaxation of the Inflaton: Consider the following Lagrangian for the inflation 
$ a scalar a and a fermion field ip: 

C=l (d<S>) 2 + ^(da) 2 -V(®,o-) + tlj(i0-M^-y<S>)tlj 

V^a) = \Ml^ + ^^ + l -Mla 2 + ^^ + 9 -a^\ (1) 

where a, ip are the lighter degrees of freedom (M$ >> M a ^M^). As in our previous work 
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H we write <&(x,t) = 5(t) + i) where 5(t) is the expectation value of the zero mode 
of $ in the non-equilibrium density matrix and (x(x,t)) = ||. The expectation value of 
the a field is assumed zero; this is consistent with the equations of motion. Following the 
formalism of 0, we arrive at the linearized, non-equilibrium equation of motion for S(t): 

6(t) + M 2 5(t) + f 5(f) E r (t - f) dt' = 0. (2) 

The kernel £ r (t — t') is obtained in the loop expansion and is determined by the retarded 
self-energy. This self-energy needs appropriate substractions which are absorbed into mass 
and wave function renormalization. From now on, all masses and couplings will be the 
renormalized ones. We choose to start the evolution at t = using the initial conditions 
5(0) = So, 5(0) = 0. This equation can be solved via the Laplace transform |j]|§. The 
inverse Laplace transform is performed along the Bromwich contour and requires the iden- 
tification of the singularities of the inverse propagator: G(s) = [s 2 + M| + E r (s)] 1 where 
s is the Laplace transform variable. We have to distinguish different cases: i) Unbroken 
symmetry case: a) If y — 0, the lowest order contribution to E r (t — t') is given by the 
two-loop retarded self-energy which has a cut for s 2 > — (M$ + 2M a ) 2 . 

For weak coupling G(s) has a one-particle pole below the three particle threshold. The 
contribution of this pole gives oscillatory behavior of 5(t) with the period determined by 
the position of the pole. The contribution from the continuum to the long time dynamics is 
determined by the behavior of the imaginary part of E r near threshold and it yields a power 
law t~ 3 at long times for all non-zero masses. The imaginary part of the self-energy does 
not indicate exponential relaxation and although the amplitude of the oscillation damps 
out, it reaches a non-zero asymptotic limit at long times. It is straightforward to prove 
that the ratio of the amplitudes at two different times is a renormalization group invariant 
||. b) If y ^ then the kernel receives a contribution at one loop from the fermion- 
antifermion intermediate state. The self-energy has a cut beginning at s 2 = — 4M|. If 
M$ > 2M.0, the one-particle pole moves off the physical sheet into the second sheet and 
for weak coupling, the spectral density has the Breit-Wigner form with a narrow resonance 
at the renormalized value of the scalar mass and width T « y 2 M^/8ir (for M$ 3> M^). 
The expectation value 5(t) now relaxes exponentially with a relaxation time t re \ = for 
t < r _1 ln(M$/r), but eventually the long time behavior is dominated by the behavior of 
the discontinuity of the imaginary part of the self energy near threshold leading to a t~ 5 ^ 2 
power law behavior. Two subtractions of the self-energy are needed, resulting in mass and 
wave function renormalization. If M$ < 2M^ there is a one-particle pole below the two- 
particle threshold i.e. in the physical sheet, and the long time behavior of 5(t) is oscillatory 
with nonzero asymptotic amplitude. The continuum contribution falls off with the power 
law above. Clearly an imaginary part for the self energy does not translate to exponential 
relaxation, ii) Spontaneously broken symmetry case: If M| = — fi 2 < 0, the new 

minimum is at $ = y6fj?J\§, and we write t) = $o + 5(t) + x(x,t). The self energy 
now has two-particle cuts beginning at s 2 = — 4(M a + g§ ) 2 , —A(M 1 p + y^> ) 2 , — 8/i 2 . If 2/i 2 is 
greater than the lowest two particle threshold the one-particle pole moves into the second, 
unphysical, sheet and the spectral density features a narrow resonance, for weak coupling, at 
the renormalized value of the scalar mass of total width given by T to t — T^, 
8n + 3^ 2 /x/(47Ta/2A$). In this case 5(t) exhibits exponential relaxation on a time scale that 
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is typically as the one above. Eventually, however, a power law tail behaving as t~ 3 / 2 , which 
is determined by the two-particle scalar threshold, will develop. 

From our discussion above, we find two important points: (i) In the case when the inflaton 
mass is above the multiparticle thresholds, the narrow resonance gives rise to exponential 
relaxation for a long period of time. During this time, the evolution is similar to that 
obtained from the effective evolution equation 5 + m 2 5 + Y tot 5 = with T tot given by the 
total decay rate. The equation of motion (0) provides a "proof" of this effective description; 
however this only applies for small amplitudes, (ii) Usually the imaginary part of the self 
energy is identified with a "thermalization" or relaxation rate. We have seen above that 
this is only the case whenever there is a resonance i.e. a pole in the unphysical sheet. If 
there is a one particle pole below the multiparticle threshold, then the long time behavior 
is undamped and oscillatory. The continuum contribution falls off with a power law which 
is determined completely by the behavior of the spectral density at threshold. 
Non-Linear Relaxation of the Inflaton: We now contrast the linear relaxation scenario 
with the non-linear one, which is found by considering evolution equation for the expectation 
value of the scalar field including the back-reaction of the quantum fluctuations by absorbing 
5(t) in the definition of the now time dependent mass term for the scalars and fermions: 
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mi(t) = M' 2 + ^ S 2 (t) , ml(t) = M 2 + g 5 2 (t) , m^t) = + y 6{t). 



(3) 



We can then obtain the one-loop evolution equations 



A 



5 + MIS + ^ 5 3 + 5 ^ ( X 2 (x , t)> + g 5 (a 2 (x, t)) + y (V(af, t)^(x, t)> = 0, 



A 
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(4) 



where the non-equilibrium expectation values are the one-loop Feynman graphs but with 
the time dependent masses given by eq.@: 



(X 2 (x,t)) 



d 3 k \U 2 {t)\ 
(2vr) 3 2^(A;) 

d 3 k 



, (o- 2 {x,t)) 



d 3 k \V 2 (t)\ 



(27T) 



1 - 2k 2 fl(t)f k (t) 



(5) 
(6) 



The mode functions Uk(t), 14 (t), fk(t) obey 
d 2 



+ k 2 + m%{t) 



dt 2 

^ + ^ + ml(t) 
d 2 

— + k 2 + mj(t) - im^(t) 



U k (t) = , C4(0) = 1 , f4(0) = 
V k (t) = , V k (0) = 1 , 14(0) = -iu a (k) 
fk{t) = o , f k (0) 



y 1 uj^k^uj^k) + m^(0)) 



/ fc (0) = -iaty(A;)/ fc (0) , uo a (k) = Jk 2 + m 2 a (0) , a = $,^,cr. 



(7) 
(8) 
(9) 
(10) 



The mode functions fk{t) define the spinor solutions of the homogeneous Dirac equation. 
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From eqs.(|7|- ||), induced amplification of quantum fluctuations becomes evident through 
the time dependent mass; the mode functions obey harmonic oscillator equations with a time- 
dependent frequency. If the time dependence were truly periodic this would be parametric 
amplification through the presence of resonances and forbidden bands |4],[5|||]. However, 
when back reaction is incorporated as in the set of equations above, these time dependent 
frequencies are not exactly periodic precisely because of the damping effects. This case 
should be distinguished from parametric amplification and we refer to this case as induced 
amplification. 

The first order quantum correction to the evolution equation is obtained by replacing the 
classical trajectories for 8(t) in the mode equations (0-0). However, this will lead to infinite 
particle production via parametric amplification. In reference |j] , the classical evolution for 
8(t) was kept in the mode equations. We, however, will keep the full evolution for 8(t) thus 
obtaining a self- consistent one-loop approximation that incorporates the back-reaction of 
the field into the fluctuations and vice- versa. In effect, this procedure is a non-perturbative 
resummation of a large set of diagrams self-consistently. 

Induced amplification will give rise to damped evolution for 8(t), whose amplitude dimin- 
ishes at long times. Thus, incorporating the full evolution of 8(t) in the mode equations will 
shut-off particle production at long times as the energy of the scalar field is transferred to 
the produced particles. This closed system of consistent equations will ensure that particle 
production will be finite. 

To obtain the dynamics of the zero mode, we will integrate the equations above numeri- 
cally, using the initial conditions 8(0) = 8q, 8(0) = 0. These conditions imply that the modes 
written above correspond to positive frequency modes for t < 0. We have to renormalize the 
equations to one-loop order absorbing the relevant divergences into mass, wave-function and 
coupling constant renormalizations || . The final form for the zero mode equation is as given 
by eq.(||) in terms of renormalized couplings and masses and the expectation values have 
been twice subtracted with the aid of a high momentum cutoff A (typically, A ~ 100|M$|). 
Such a cutoff is necessary in any case in order to perform the momentum integrals appearing 
in the expectation values numerically. The results have been checked to insure that there is 
no cutoff sensitivity in physical quantities. The higher mode equations take the same form 
as in eqs.(^|- |]), but in term of the renormalized parameters ||||. 

We now use these equations, written in terms of the dimensionless quantities r](r) = 

$(t)^/A$/6M|; r = M$ t, to anyalze the dynamics of the zero mode as well as to com- 
pute the number of a, ip quanta produced by induced amplification ||6|||. iii) Unbroken 
symmetry case: a) y — 0: Figure 1 shows rj(r) vs. r for A$/87r 2 = 0.2, g = A$, M a = 
0.2 M$, 77(0) = 1.0, 77(0) = 0. The amplitude damps out dramatically within a few periods 
of oscillation. The asymptotic behavior corresponds to small undamped oscillations. This 
is to be contrasted with the linear relaxation case, in which the two-loop self-energy leads 
to undamped asymptotic oscillations and a power law behavior at long times. We also find 
that the number of particles in a volume M$ 3 , as a function of wave number k, is peaked at 
k < M$. Induced amplification is much more effective than one-particle decay in producing 
low momentum particles and drawing energy from the background field, thus damping its 
evolution. One-particle decay considerations would predict a slow power law tails to a final 
amplitude that is perturbatively close (notice that the perturbative parameter is \<p/8ir 2 ) 
to the initial. 
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We also find || that both the quantum fluctuation (a 2 (x,t)) and the total number of 
produced particles as a function of time grow rapidly within the time scale in which damping 
is most dramatic, as can be seen in figure (1), clearly showing that damping is a consequence 
of induced amplification of fluctuations and particle production, b) y 7^ 0, g = 0: Figure 
2 shows rj(r) vs r for M$ = 0, y 2 /^ 2 = 0.5, A$/6y 2 = 1.0. Clearly there is some initial 
damping that shuts off early on. This is in marked contrast with the prediction of the linear 
relaxation approximation that predicts exponential relaxation with a relaxation time (in the 
above units) r re i « 0.1 for the above value of the couplings. We find || that the number 
of particles in a volume M$ 3 as a function of wave number is peaked at k < M$ with 
maximum value of 2 given by the exclusion principle. We interpret the lack of damping as a 
consequence of the Pauli blocking; all the low momentum states are fully occupied. In the 
initial stages of evolution an excited state quickly ensues and the available fermionic states 
are occupied blocking any further damping. Both the quantum fluctuation (^)(x,t)ij;(x,t)) 
and the number of produced particles as a function of time grow very rapidly within the time 
scale of the initial damping, but then oscillate around a constant value at the time when the 
damping effects have shut off ||. iv) Spontaneously broken symmetry case: y = 0. 
Figure 3 shows the time evolution of t)(t) vs t for A$/87r 2 = 0.2, g/\<s> = 0.05, M a = 
0.2 M$, 7/(0) = 0.6, 7y(0) = 0. For these values of the parameters, linear relaxation predicts 
exponential damping with a relaxation time r re i ~ 360. This relaxation time is about six 
to eight times larger than that revealed in figure (3). Furthermore, this figure clearly shows 
that the relaxation is not exponential. The asymptotic value of 7/(r) oscillates around the 
minimum of the effective action which, as observed in the figure, has been displaced by a 
large amount from the tree level value rj = 1. We notice numerically a slow damping in this 
asymptotic behavior, which could be related to the linear relaxation time, though we were 
unable to establish this numerically. We find [^] that the number of produced particles in a 
volume M$ 3 is peaked at low wave number with a fairly large contribution within the region 
k M$. As a function of time the quantum fluctuation (o~ 2 (x,t)) grows rapidly within 
the time scale in which damping is most pronounced and the same behavior is observed for 
the number of particles created as a function of time, thus leading to the conclusion that 
the rapid non-linear relaxation is clearly a result of particle production and the growth of 
quantum fluctuations. In all cases we find the total density of produced particles to be 
M = AMf , JVwO(l). 

Conclusions: We have learned that induced amplification of quantum fluctuations in the 
non-linear regime is a highly efficient mechanism for damping the evolution of the inflaton 
zero mode as well as for creating particles. The typical time scales associated with this 
damping mechanism are much shorter than the time scales associated with single particle 
decays. This leads us to conclude that there may be a wide separation between the time 
scales of induced particle production and that of thermalization of these produced particles 
via collisional relaxation, since this latter scale is typically related to perturbative couplings 
and thus fairly large. Single particle decay will be operative on the same time scale as 
collisional thermalization. 

The final reheating temperature, however, will be determined by the total number of par- 
ticles produced. Typically, T re h ~ A/" 1 / 3 . Similarly, the total entropy produced will depend 
on the number of produced particles. Assuming that the process of single particle decay 
produces far fewer particles than induced amplification and that a thermalized state ensues 
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from collisional relaxation, one can estimate the final reheating temperature by evolving the 
equations above and following the evolution of the number of particles. Clearly this number 
will depend on the particular details such as coupling, masses but not in the perturbatively 
obvious manner. Assuming that the produced particles thermalize via collisions in the cases 
that we studied we typically obtain reheating temperatures T re h ~ M$ [§] which is about 
two orders of magnitude less than the estimate from the elementary theory of reheating [|TJ . 
However, it is not clear how long a distribution that is so different from thermal will take to 
thermalize. In an expanding universe, this time lag will serve to dilute the energy density of 
the produced particles. Ultimately, a realistic calculation will involve the numerical evolu- 
tion of the self-consistent equations. Furthermore our study will have to be extended to an 
FRW cosmology since the expansion of the universe can lead to many effects not accounted 
for here. This work is in progress. 
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Figure Captions 

Fig. 1: Unbroken symmetry case: y — 0, A$/87r 2 = 0.2, g = A$, M a = 0.2M$, 
1.0, t)(0) = 0. 

Fig. 2: Unbroken symmetry case: g = 0, = 0, y 2 /V 2 = 0.5, A$/6y 2 = 1.0, 
1.0, t)(0) = 0. 

Fig. 3: Broken symmetry case: y — 0, A$/87r 2 = 0.2, g>/A$ = 0.05, 
0.2 M$, 77(0) =0.6, t)(0) = 0. 
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